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Abstract 



A theory of superexchange in the mixed valent layer compound NaV205 is 
presented which provides a consistent description of exchange both in the dis- 
ordered and charge ordered state. Starting from results of band structure 
calculations for NaV205 first an underlying electronic model for a ladder unit 
in the Trellis lattice is formulated. By using the molecular orbital representa- 
tion for intra-rung electronic states a second-order perturbation procedure is 
developed and an effective spin-chain model for a ladder is derived. Variation 
of the resulting superexchange integral J is examined numerically as the lad- 
der system evolves from a charge disordered to the extreme ('zig-zag') charge 
ordered state. It is found that the effective intra- ladder superexchange is 
always antiferromagnetic. 
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I. INTRODUCTION 



Since the phase transition into a spin-gapped phase in a'-NaV205 was reportedB this 
layered insulating compound has attracted much attention. It was found that the phase 
transition (T^ ^ 34 K) is accompanied both by the lattice distortion! and a charge ordering! 
in the vanadium layers. Though the details of the low-temperature phase structure in 
NaV205 remain controversial, some main features of this structure are obtained from ex- 
periment. For instance, the same low-T lattice suDerstructure was measured and reported 
by three groups of researchersSH^. The experimentQH has also ruled out a hypothetical 'in- 
line' charge ordering (CO) in favour of a 'zig-zag' CO. At the same time mechanism of the 
spin-gap formation is not understood up to now and formulation of adequate spin models 
for NaV205 is required. 

The electronic and magnetic properties of NaV205 are mainly due to the V ions having in 
average a mixed valence -1-4.5. In a particular layer V-ions are arranged to give a sequence 
of two-leg quater-filled ladders (one electron per rung) coupled to each other in a Trellis 
lattice. It is expected that the main properties of vanadium layers are dominated by the 
ladder physics and the theoretical consideration should start with an analysis of this one- 
dimensional structure unit. The most important evidence in favour of this point of view is 
provided by the band structure calculations^^ within the density-functional theory (DFT). 
These calculations have shown that the overall electronic band structure of NaV205 is 
determined mainly by the large intra-ladder electronic hopping amplitudes while the inter- 
ladder hopping is rather weak. 

Let us now turn to the problem of spin coupling in NaV205 and consider a separate 
quarter- filled ladder. Since the strongest electronic hopping transfers ~ t^ between the 



basic vanadium d^^^ -orbitals are within the rungsaliy one has to ascribe a spin-1/2 not to a 
vanadium site but to the intra-rung V-O-V bonding molecular orbital. With strong Coulomb 
repulsion ~ f/ on the dxy -orbital the lowest electronic bonding orbital cannot be doubly 
occupied and the problem is naturally reduced to a half-filled chain in the limit of strong 
electron correlations. The excitation energy of the lowest 2-particle singlet rung state from 
now on plays a role of an effective Hubbard repulsion parameter which is roughly estimated 
to be 2t(j in the charge disordered statei0. Horsch and Mack were the first0 who derived 
a quantum Heisenberg spin-1/2 chain model to describe a one dimensional behavior of the 
ladder systems NaV205 in the charge disordered high-T phase. 

The superexchange coupling parameters entering into the effective Heisenberg chain 
model are expectecOo to be strongly reduced as the system undergoes the phase transition 
into the low-T charge ordered state. It was predicted, for instance, that for neighbouring 
spins in a ladder the exchange integral J, which is antiferromagnetic in high-T phase, may 
even change the sign to become a ferromagnetic one in the extreme CO phase with a zig-zag 
charge distribution on the ladder rungsl^. In the present paper we develop the superex- 
change theory and obtain estimates for the spin-spin coupling constants in a single ladder 
both in the charge disordered and in an ordered phase. That allows us to reexamine some 
earlier statements and estimates0^El. Since the in-line CO is ruled out by experiment the 
present consideration is restricted to a zig-zag CO state only. 

The nature of the CO phase transition in the real NaV205 compound is not well under- 
stood up to now. Nevertheless, at the first stage (Section HI) the analytical expressions for 
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the superexchange integrals in this system can be derived and analyzed in a rather general 
context. Next, to reach a quantitative description one has to rely on a particular realistic 
mechanism driving the charge ordering. We assume that the intersite Coulomb electron 
repulsion is responsible for the CO phase transition and the lattice degrees of freedom play 
a secondary role. The inter-ladder short-range Coulomb interaction is treated in a mean 
field approximation while intra-ladder charge-charge correlations are considered on a more 
accurate level. To this end considering at some stage (Section IV) the charge degrees of 
freedom only we transform the electronic Hamiltonian onto the effective Ising model in the 
transverse magnetic field0@. The pseudospin variables in this approach are associated with 
the charge dynamics and the order- disorder transition is determined by the dimensionless 
interaction parameter A = Vb/2ta- Here Vb is the characteristic Coulomb repulsion be- 
tween electrons on neighbouring rungs. We adopt this approach which allows us to examine 
quantitatively how the superexchange coupling varies as the ladder system evolves from a 
disordered to the charge ordered phase. We expect that besides the quantitative estimates 
this particular analysis reveals general, mechanism-independent, trends in the behaviour of 
the superexchange spin-spin coupling in NaV205 . 



II. ELECTRONIC MODEL AND THE EIGENSTATE PROBLEM 

We start with the electronic Hamiltomian for a quater-filled ladder in the following 
general form 

H = Y. {h^' + + E {H^f + H^'') (1) 

i <ij> 

where (i) and (ij) refer to a particular rung or to a pair of nearest neighbouring (nn) 
rungs, respectively. The intra-rung interactions involved into {H^^ + H^^]_) are described 
by the electronic hopping amplitude ta, the Coulomb on-site, U, and inter-site, Va, param- 
eters. For the nn inter-rung coupling {H^^^ + Hy'''^), there are two, tb and t^, electronic 
hopping amplitudes and the corresponding Coulomb repulsion is given by Vb (see Fig.l). By 
introducing the creation operator dlj^^ {dln^) ^ spin- a electron on the left (right) V-ion 
within the i-th rung we specify the particular terms in Eq.(|I]) as follows: 

H^'^ =eY,K + V^n^nf + U <r<l (2) 

a=L,R a=L,R 

-ta (dlLadiR^ + h.C.^ (3) 

cr 

E E {4aAa. + h.C.') (4) 
a=L,B. a- 

-td {dlladjRa + 4j^djLa + k.C.^ 

(T 



a=L,R 
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where nf = n"^ = J2a dlaa^iaa (tt = L, R) and e is a bare vanadium d^^y-orbital energy 
in an isolated ladder. 

Let us fix now the range of the model parameters. We use the following set of hopping 
amplitudesS: ^^=0.38 eV, th ~ td =0.085 eV. The proper signs of these amplitudes are 
adopted in the definition of H. Importance of additional hopping processes ~ td, which 
couple dxy- orbitals of two V- ions at the opposite ends of nn rungs of a ladder (Fig. 1), was 
reported in Ref.( [1^). Approximate equality td — tt, allows one to explain the pronounced 
feature of the band structure in NaV205 . Namely, there is almost no dispersion of the 
antibonding- type molecular orbital states formed by the d^y- orbitals on a rung. We will 
consider also another set: ta=0.3S eV, tb=0.17 eV and td=0, as a more representative one 
among the sets with td=0 used by other authorsiJnFEl. The value of the on- site Coulomb 
repulsion U is fixed from the supercell LDA+U calculations^ to be f/ = 4 eV. The estimates 
for the magnitudes of intersite Coulomb repulsion Va and H (Fig. 1) given in the literature 
are less accurate. In Ref.( |13]) for instance, it is claimed that a screening reduces Va and 
Vb parameters to rather small values Va, Vh < 0.4 eV, and hence the inter- rung Coulomb 
repulsion cannot stabilize CO. To investigate this point we developed a combined analysis 
by using the results of DFT calculation^ and more generic multi-orbital Hubbard model for 
a ladder in the Trelllis lattice. We found that inter-orbital vanadium d-d electron correla- 
tions lead to an enhancement for the effective V^ parameter which enters into the effective 
single-orbital Hamiltonian Eqs.(|I]-^. These results will be presented elsewhere^. In the 
present investigation we take an intermediate value for Va —0.5 meV and allow Vb to vary 
in a wide range up to 1 eV. 

It is worth noting that the validity of the perturbation procedure developed below to the 
second order in the inter-rung hoppings is based mainly on the smallness of these hoppings, 
tb and td, compared to 2ta provided the on-site Coulomb repulsion U is the largest param- 
eter of the theory. This guarantees also that the quarter- filled ladder is in an insulating 
stateB0lli. In contrast to the perturbation approach developed in Ref.( |1^) in our case there 
are no formal restrictions on the values of Va and Vb- Like in Ref.( 0) these parameters 
could be arbitrary small giving a meaningfull description of the superexchange in a charge 
disordered state of NaV205 . Surely, Vb should be large enough, Vb > V^^^^\ to sustain a 
zig-zag CO state and the critical value V^ depends on the mechanism chosen for the CO 
phase transition. 

In the perturbation procedure the zeroth order Hamiltonian is defined a.s Hq = H — 
'^^ijy "^ty "^here the last term is a perturbation mixing electronic states on nn rungs. As 
the necessary prerequisite of the superexchange theory the 1- and 2- particle electronic states 
on a particular rung have to be found explicitly. That can be done by considering the inter- 
rung Coulomb interaction, the term Yl<ij> Hy'^^ in Hq, in a mean field (MF) approximation. 
In this case the inter- rung Coulomb interactions lead to the following shift of the on- site 
orbital energies (a = L,R): 

e.o. = e + Vb + {nU)g\ (6) 

where {nfj^^)g are the averages over the lowest singly occupied electronic states of two rungs 
neighboring to the i-th rung within the same ladder. Having ^^('^f)^ = 1, one obtains 
{SiL + £^iK)/2 =e + Vb = e and 
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-1)*A£; Ae = 2Vbr] 



(7) 



In Eq.(0) we introduced the parameter rj {0 < rj < 1) oi charge disproportionation 
{niR)g — {niL)g= {—^Yv ill zig-zag CO state of a ladder. 

After applying the approximation given in Eq.(^ the zeroth order Hamiltonian Hq be- 
comes a sum over individual rungs, Hq ~ ^i-^^o*''- Fro^i now on each term Hq^ can be 
diagonalised independently. It is worth emphasizing that in this approach the local intra- 
rung correlations due to U and Va are treated exactly. To specify explicitly a form of Hq^ 
in a particular n— particle sector of the intra- rung electronic states we invoke an extra sub- 
script like in HqI^ below. For instance, by solving the eigenstate problem in the 1- particle 
sector one obtains 



where the energy Eg/f of the lowest /excited molecular orbital is given by 

eg/f = ETl/2ViAsr + {2tay 

and 



(8) 



(9) 



u. 



(10) 



The Ui and Vi coefficients in Eq.([TO|) are defined in the following way 



\ 



1 + 



Ae,; 



^{Aey + (24) 



\ 



Ae, 



^{Aey + {2ta) 



(11) 



and for Asi ^ the values of u, and Vi alternate as i runs along the ladder. The complete 
basis for the 2- particle states on the i-th rung consists of three singlets \Si^n > {n = 1, 2, 3) 



and triplet states Ir/™'' > (m = 0,±1). Explicitly, we choose the following definitions for 
these states (the rung index i is implied): 

l5i) = i/v^(44i-44T)io) 
1^2,3) = i/v^(±4A + 4T4i) |o) 

|rW) = l/v^(44^ + 4,4,) |0) 
|r(±i)) = 44|0), 4i4i|0) 



(12) 



One can check that the triplets are degenerate eigenstates of -^0*25 i-^- Hq2\t^'^') = Erlr-"^'), 
with the energy er= 2e + Va. In the singlet subspace the Hamiltonian H^l takes the following 
form 



r(«) uM\ 
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(2e + Va -2ta 

Hl,i = {\Si,i),\Si,2)AS,,3),)\ -2ta 2e + U Asi | | (5^,21, | (13) 

\0 Ae, 2e + U 

By diagonalizing the symmetric 3x3 matrix one obtains the eigenstates 

l^f^)=E7i^n^)|5.n) (14) 
n=l 

and the corresponding eigenvalues Esk ('^=1,2,3) in the singlet subspace. More precisely, 
= 2£ + Va+ {U — Va)eK, wherc e^in = 1, 2, 3), are the solutions of a cubic characteristic 
equation reading 

e4{l-e^y-C^]+l2^[l-eA = (15) 

where Q = Aei/{U — K) and /i = 2ta/{U — Va). Note, that according to Eq.(0) the value 
of C,f does not depend on i. Provided Eq.([l5|) is solved the 3-component vector in 
Eq.(|np can be found by using the following expressions 



7i 



lA 



7r) = ±.z^n(A.)^^ii-^ (16) 

where = {1 - e^f - Q and R,, = [Al + /^^(l - e^f + jj'^CfY^'^. Note, that the upper/or 
the lower signs in Eq . ([T7|) should be taken simultaneously. Like Q, the third component 
73 (i) is staggered along the chain direction. This particular problem, Eqs.(|15|, [1^), can be 
solved numerically in a wide range of the underlying model parameters provided the order 
parameter t] is fixed from some self-consistent procedure. Most representative results will 
be discussed in Sec.V. 



For further purposes it is helpful to use the projection operators Xf''^ = \qi){qi\ and the 
transfer operators X/'^ = \qi){pi\ where \qi) and \pi) are the states defined above on the i- 
th rung. More explicitely, for instance, one has 

Xr'' = \g.„){g.^\, X^'''' = \0,){g,^\, (17) 
^Sk,Sk ^ .9<x ^ \sl^)){g^^\ etc. 

With these notations the approximate zeroth order Hamiltonian Hq ~ - Hq^ can be written 

as 

H^^ = e.Xr + e,Y. Xr" + 6/ X/-^'^ + 

(J (7 

3 

K=l m=±l,0 
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In Eq.([T8|) the vacuum rung state |0j) with the reference energy is involved for complete- 
ness. 

At the same time the 3- and 4- particle states are dropped as these states are irrelevant 
for the present purposes. At the end of this section the following remark is worth to be 
made. Actually, while solving the intra- rung eigenstate problems we could not avoid the 
use of a MF treatment of inter- rung Coulomb interactions, the term Y2<ij> Hy''-' in the 
original Hamiltonian of Eq. . Nevertheless it does not mean the complete neglect of the 
many body effects due to Hy''^ in our approach. Further analysis in Sees. IV,V will show 
that charge-charge correlations between singly occupied nn- rungs are important to obtain 
the estimates for the superexchange coupling between spins in a ladder. 



III. PERTURBATION PROCEDURE AND 
SUPEREXCHANGE ON A LADDER 

The superexchange coupling between two singly occupied rungs is established due to the 
electron hopping transitions with the virtual intermediate 2-particle states in one of two 
rungs. These transitions are due to the kinetic part of the Hamiltonian, y, whose action 
has now to be represented in terms of the molecular orbital states \qi) defined in the previous 
section. In a rather general form one may write: 

= E = E E E i\p'.p)xf''xf'^^ + h.c] (19) 



<ij> <ij>,a qq' pp' 



where 



tijA(i'^(i\p'^p) = -tb'^iq'M^Mip'jldjaalPj) 

a 

-td [{q[\d\n^qi){p'M3L^P3) + (R ^ ^)] (20) 

and {R ^ L) means the indices interchange. All the necessary amplitudes {q[\d\^^\qi) and 
their conjugates can be now straightforwardly evaluated to give the relevant part of Hi^\\ in 
the following form 

= - E ^4 ^ E (xf-^x;-- + h.c) + 

ij,a K. K=l 

+E E B^^c^jrj [xr^-' x^^^^ + h.c.)] (21) 

cr' m=±l,0 J 



Here ^ij = 1 if (ij) is a nn pair and equals zero otherwise; {2a) = +l(a =t), — l(a =|) and 

iefi 

Cirj = 5.',.5,n,2<x + l/^/26„',,6m,o. (22) 



a = —a. The factor C^™"* is defined as 



The effective transfer amplitudes and Bij can be presented in the following form: 
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n=l 



Bij — thbij + tdbij (23) 
where ajj,„, ajj^„, bij and b^ are bilinear functions of the u-, v- coefficients 

aij,i = aij^2 = UiVj + ViUj 

aij,2 = Oij, 1 = UiUj + ViVj (24) 

«ii,3 = hj = -UiUj + ViVj 

bij = aij^s = -UiVj + ViUj 

To derive the superexchange spin couphng to the second order in the effective hopping 
amphtudes d\'^'' and Bij we use the Schrieffer- Wolff transformation H'=exp(— S)-?/ exp(S') 
with the generator S{= —S'^) being determined from the condition [Hq,S] = —Ht^\\. Then 
the second order correction to Hq is given by 

Feff = -l/2[5,i/i,||] (25) 

From Eq. (^) we obtain first the form of the generator 

3 r)M 




h.c. ) (26) 



One can see, that there are two perturbation channels, (s) and (r), with the singlet and 
triplet intermediate virtual states, respectively. Finally we obtain from Eq. (^51) the superex- 
change Hamiltonian 

<ij> 

with the exchange constants in the singlet and triplet channel given by 

j{r) ^ B^ + Bl ^^^^ 

In Eq. (p^) the spin- | operators 5"^ = 1/2^^^, \gia)o'aa'{9io-'\ and a"^, are the Pauh matrices 
{a = X, y, z); the notation < ij > means that each pair of the nn rungs is taken only once. 
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Note that already at this stage the expressions Eq.(]28|) and Eq.(^) together with the 
definitions in Eqs. (0), (p4D, formulate a complete scheme to estimate j'^^^-^ and j''^^^^ both 
in a disordered and a zig-zag CO state. That can be done, at least numerically, by solving 
first the eigenstate problems of Sec. II for given parameters of the model and with the order 
parameter rj varying in the range < r/ < 1. In this approach the effects of charge ordering 
upon spin coupling in a ladder can be examined in its main features without specifying a 
mechanism driving the CO transition in a quarter- filled ladder system. Preliminary analysis 



shows that the total superexchange coupling J(jj) = J^^^-^ + J^^^ is reduced as the ladder 
develops from a disorderd to a zig-zag CO state. That is in accordance with a predicition 
in Ref.( |T2p. A more satisfactory quantitative description could be done if one relies on a 
specific mechanism responsible for CO phase transition. In the next section we pursue one 
possibility by noting that the underlying model, Eqs.(|I]-|^) contains the expected instability 
without invoking extra degrees of freedom like charge-phonon couDling. Further analysis is 
based on the pseudospin description of the charge correlationsljil^. Before turning to the 
pseudospin model itself (Sec. IV) we use this representation to present the main results of 
this section given in Eqs.(|2^, 0, |2^, in a more condensed form. 

In the pseudospin representation the molecular ground state of the i-th rung \gic) = 
\9i) ® is given by 



2 I 2 



(30) 



By using the pseudo- spin (T = |) operators T^'^a = x,y,z) we define now several pair 
correlation functions and expectation values in the ground state of a ladder as follows 



Q 
R 
R 



(x) 
{xz) 



{T^ + T,-),, R 



(^) 
j 

(^) 



(31) 



By using the above definitions and Eqs. 
expressions of Eqs.(E3, E9]) in the following form: 



53) one may write the numerators in the 



D 



+ 7i 



Bfj + (^ - j] 



7r«-7r(j)|/Cr+7^' 



. n=l 
^(13) 



2r-(4) 



7r^(0 



(32) 
(33) 



where 



1^ 
IC 
/C 



(1),(2) 
ij 

(3),(4) 
ij 

(12) ^ 
ij 



2v 



(1 + .^) 
(1 + v^) 

1 



1 
4 
1 
4 

(^) 



T (1 - y') + 2uR 

± (1 - of 

+ (1 + R 



(x) 



(34) 
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^(13) ^ _ 



and we defined u = t^/tf,. It is worth noting that the antisymmetric property, /C^j^'' = — /C 
and /Cj-^^''= — /Cj^'^\ makes the complete expression of Eq.(^) to be symmetric. 

IV. PSEUDO- SPIN MODEL FOR THE LADDER 



(13) 
'ji 



The effective Hamiltonian in Eg . (1271) describes spin-spin interactions in the low- energy 
subspace spanned by the vector charge- spin manifold |G'{g}_{o-})= HJIS'*) ® l^*)]- ^.im 
now is to project the electronic model, Eqs.(|l]-^ onto the charge sector within this manifold. 
The projection procedure was discussed several timeJiiJii, it results in the following zeroth- 
order effective Hamiltonian for a ladder 



= 2H [TIT! + 1/4] - 

i <ij> 

-Ae'J2^-l)X' (35) 



In this one- dimensional (ID) model the effective longitudinal field ~ expresses in a 
condensed form effects of Coulomb interactions due to higher dimensionality of the real 
NaV205 compound (see Appendix). 

We proceed further in the spirit of the coupled quantum spin chain approach0 which 
involves a MP treatment of the inter-chain coupling and a more sophisticated treatment of 
the intra-chain, i.e. intra-ladder, interactions. The intra- ladder problem enters into the 



expressions for the superexchange constants, Eqs.(p8|, 29) and Eqs.(31-p^, in the form of 
pair correlation functions for the nn pseudospins. Having the final aim of calculating these 
constants we rely on properties of the exact solution for the ID Ising model in a transverse 
field. 

The term ~ in Eq.p5|) breaks explicitly the Z(2) symmetry of the ID transverse 
Ising model which results in a non-zero longitudinal magnetisation {T^)g = {—iyr]/2 if the 
coupling constant A = Vb/2ta exceeds some critical value, A > Ac. One may write 

(36) 



A 

and choose for the exponent (3=1/8 since the presence of a weak longitudinal field does 
not change the universality class of the ID Ising model in the transverse field. Note, that 
for this exactly solvable model Ac = 1 (Ref. 0,0) and we keep this value in the subsequent 
numerical analysis. The open problem that still remains are unknown dependencies of f3 
and Ac on details of the inter- ladder Coulomb interaction if one goes beyond a mean field 
treatment of these interaction in the Trellis lattice. 
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Let us now turn to a discussion of the superexchange constants, Eqs. (p8|, pQl). We estimate 



the inter- rung pseudospin correlation functions entering into Eq.(31) at two different levels 



of sophistication. The simplest one is to decouple the pair correlation functions as follows: 

Qlf^(^^)' = 1/4(1-^^) (37) 



In Eq.(|37|) we take into account that i and j are refered to nn rungs and, hence, (Tj) 



— {Tf) while {Tj') = a/1 — ^7^/2 does not depend ont the rung index i. Within this first 
approximation the superexchange constants will be calculated in the next section. The only 
signatures of the exact solution in Eqs.(^, are due to the choice of the values for 
and Ac, while the inter- rung fluctuations are ignored. To gain an impression how these 
fluctuations may change the superexchange constants we adopt some more properties of the 
exact solution for the ID transverse Ising model. For this purpose let us considerEl'El the 
pair correlation functions {j = i ± 1): 

{TIT^), = K- {T:T^), = ml + K- {T^Ty), = ky (38) 

In particular, one haSlli 



where Ag = y 1 + + 2Acosg. The expressions for k^^y can be found in Refs.( p!8| , |l9D as 
well. We use the property jfej;!, -C rrix, \ kz\ and drop the weak transverse fluctuations, 
which means k^^y — 0. Breaking of the Z(2) symmetry allows us to write also 

^(-) ^ (40) 

We use the Eqs. (|36|, pH]) and (|39D together with the approximation in Eq.(|40|) as the basis 
of our second approach to calculate the superexchange constants. 



V. ANALYSIS OF SUPEREXCHANGE INTEGRALS, 
NUMERICAL RESULTS 

The starting point of our analysis is the exact diagonalisation of the single-rung Hamil- 
tonian in the 1- and 2-particle sectors. The excited singlet and triplet states provide the 
intermediate states for the superexchange mechanism of magnetic coupling between neigh- 
bouring singly occupied rungs. From now on we drop the irrelevant lower indices at J. The 
excitation energies Eg^ = Esk + £v — '^^g and Er = Sr + e.^ — 2eg which enter as the denom- 
inators into the perturbative expressions for J^''\ J^'^\ are determined by the eigenvalues 
due to Eqs.(P) and (^). These eigenvalues depend on the degree of charge order via the 
intra-rung d- orbital levels shifts Ae: which is proportional to the CO parameter rj in Eq.(8). 
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The symmetry breaking source is assumed to be arbitrary small. The excitation energies in 
the singlet and triplet channels as function of Ae are shown in Fig. (2), = 0. There is a 
crossing of two singlets at very high values of A^j. For realistic values of V;, however this 
crossing point is not reached. 

While solving the eigenvalue problem in the singlet sector the components of the vec- 
tors Eq.([T7|), are calculated as well. These vector components together with pseu- 
dospin correlation functions determine the numerator, Eg. (P^) , in the expression for J(^), 
Eg. (pHD . The pseudospin pair correlation functions are calculated either in the modified 
MF approximation, Eg. (P7|) , or by using a more accurate approach based on Egs.(^Sp-(^Dp. 
In the latter case the most pronounced charge fluctuation contributions to the superex- 
change are involved. This reguires evaluation of only two pair correlation functions, namely, 
{TfTf^^)g = h and {TfTf^^)g ~ ml which are given by Eg.®,®. Both of them are 
shown in Fig. (3) as function of A~^=2ta/Vfe. It is seen that the intra- chain fluctuations 
of the order parameter described by the 'longitudinal' correlation function extend far 
above the critical value A^^=l where the CO vanishes. Since in our approximation k^, ~ 
the 'transverse' correlation function is reduced to the sguare of the 'covalency' of each rung 
rrix = {Tf)g. The latter approaches one half for oo. It shows the opposite behaviour 
compared to the complementary charge order parameter = \ {Tf)g\ = r]/2. The latter is 
shown in Fig. (3) as given by Eg.(p6D with the exact ID Ising exponent (5=1/8. 
With these necessary ingredients it is now straightforward to calculate the superexchange 
integrals J^^\ J^^^ as function of the microscopic model parameters according to Egs.(p^), 
(^). The hopping and interaction parameters fixed from LDA+U- calculations are discussed 
in Sec. II where definite values for ta,tb,td and U were given. The intra- rung repulsion Va 
was assumed to lie close to 0.5 eV and the inter- rung repulsion VJ, is allowed to vary in a 
wide range up to 1 eV. The dimensionless control parameter for CO is X=Vi,/2ta with Ac=l 
at the CO transition. 

The dependence of J^^\ J^^-* on Vb is shown in Fig. (4) for two cases: The dashed line 
corresponds to exchange integrals calculated in the modified MF approach. Eg. (pT]) , with 
the Ising exponent (3=1/8. For A < 1 (H < V^*''^'^''' = 2ta = 0.76 eV) in the disordered ladder 
the exchange constants are independent of Vb and the total value J = J^*-* + J^'^^ amounts 
to Jmf = Jmf — 80 meV since Jmf^ = in this region. At Vb = v/*^^'' a sharp drop in the 
singlet exchange constant is observed whose shape is determined by the order parameter 
rriz (Fig. (3)). The most important reason for the reduction of superexchange is the reduced 
effective hopping between adjacent rung bonding states caused by charge order. The full 
line in Fig. (4) shows the same exchange integrals with the effect of charge fluctuations 
included, as described by the pseudo- spin formalism of the previous chapter, Egs.(|38D-(^0D. 
In the disordered regime the inclusion of charge fluctuations strongly reduces the exchange 
constants compared to the MF value. The total superexchange integral falls now into the 
region 70 meV> J > 55 meV, which is somewhat lower than the theoretical value obtained 
in Ref.( [Tl|) for the disodered phase while it is somewhat higher than the value j(^-'^P) ~ 50 
meV estimated in Ref.( |l]). Fig.4 shows that in the CO regime {Vb > V^*^^^), where charge 
fluctuations are strongly suppressed, both approximations give very similar results. Above 
the critical region at, for instance, Vb — 1.3V^^*"^^ (i.e. at VJ, = 1 eV), where CO is nearly 
complete, ruz ~ 0.43, the superexchange integral reaches the value J ~ 0.22 meV. 
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In Fig. (5) we show J^^^ and J^'^^as function of V^, with charge fluctuation effects included, for 
various intermediate values of the intra-rung Coulomb repulsion Va- Moreover, by varying 
Va in a wider range from zero up to rather high values, Va ~ 4ta, a regular behaviour 
of the exchange constants is found, which is already seen from Fig. (5). First, the triplet 
contribution J^^^ is almost negligible as compared to the singlet part J*-*-* for the model with 
td = tb- In the case td =0 however J*-^-* gives a sizeable (negative) contribution to the total 
J. Secondly, with increasing Va, the total superexchange integral J decreases and remains 
antiferromagnetic (> 0), it never changes its sign. This is true for both models considered 
{td = h or trf = 0). These observations are at strong variance with the results of Ref.( |13|) 

It is instructive to examine also both the limit of the extreme CO, 1/2, which 

is reached at A^^ 0, i.e. at unphysicaly high values of H, and the the opposite limit 
— > oo, i.e. at V^ = 0. To check these limits a variation of the total superexchange 
integral J is examined in a wide range of A~^, Fig. 6 (a). Two curves there correspond to 
two sets of the hopping parameters used: (1) ta = 0.38eV, U = td = 0.085eV and (2) 
ta = 0.38eV, tb = 0.17eV and td = 0. The second set is the most representative one among 
the sets with td = used by other authorsUHhEl. First, at A~^ — > the curve (1) reaches 
the value J = 7.2meV which coincide exactly with the estimate based on the standard 
expression J = 4:tyU. The latter estimate is due to the superexchange coupling between 
two vanadium ions located at the opposite ends of nn rungs. From Fig. 6 one can see that in 
the case of finite td the superexchange is regularly enhanced as compared to the case with 
td = 0. In fact we found numericaly that if A^^ oo, the two dependences, i.e. the curves 
(1) and (2) in Fig. 6, approach the same value J = SOmeV which is nothing than the MF 
value of J depicted in Fig. 4. In this limit (H = 0) the general formulas, Eqs. (|28|) , (|29|) , can 
be easily treated analytically to give the following expression for the total superexchange 
integral 

j = 2t^l(l±j£ I I 1 .41) 

I l + /i2 2ta-Ja + Va I + fl^ 2^ + J, + f/ J ^ ^ 

where /i = 2ta/iU - Va); Ja = {2taf/{U - Va) and t = t^ + td for the set (1) or t = 4 
for the set (2). By using the same parameter values as in Fig. 6 one finds that Eq.(^TD 
reproduces completely the above mentioned numerical value for J. A large change in J 
seen in Fig. 6 (a) at A = Ac is provoked by a sharp behaviour of the order parameter 
as the system undergoes the CO phase transition. Considering now the CO regime only 
(A^^ < 1) we reproduce in Fig.6(b) the same variation of J, however, with respect to the 
order parameter m^. Fig. 6(b) shows a rather smooth behaviour with almost linear decrease 
of J in the regime of weak charge disproportionation, < 0.2. 



VI. DISCUSSION AND CONCLUSION 



In Sec. II we formulated the underlying electronic model for the compound NaV205 with 
most of the model paremeters (ta, tb, td and U) fixed from known results of the DFT calcu- 
lations. By using this model in Sec. Ill we developed the superexchange theory and derived 
the second-order perturbation expression (in the effective V-V hopping amplitudes tb and 
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td) for the intra-ladder exchange integral J both in the charge disordered and in the zig-zag 
charge ordered state. Assuming that the Coulomb mechanism drives the charge ordering 
(Sec. IV) the exchange integral J is calculated in a wide range of the variable H and Va 
parameters. Since qualitative results are presented in details in Sec.V we summarize here 
most important qualitative findings. 

(i) In the intra-ladder superexchange mechanism the triplet channel is of minor impor- 
tance for the model with = and the resulting exchange integral is dominated by 
the singlet channel, J ~ J*^"*^ (Figs. 4, 5). For the model with td = 0, J^'^^ also gives a 
sizeable (negative) contribution to J. 

(ii) In the charge disordered state inter-rung charge-charge correlations within a ladder 
cause a strong reduction of the exchange integral value (Fig.4). 

(iii) With increasing of the zig-zag charge ordering the exchange integral J reduces strongly 
but never changes sign, i.e. J remains antiferromagnetic (> 0) for the whole parameter 
range. Figs. 4-6. 

(iv) Variation of the underlying electronic model due to different choices of the hopping 
amplitude (either 7^ or = 0) may lead to a considerable quantitative change 
in the value of J, especially in the extreme CO limit where J stays finite or vanishes 
respectively (Fig. 6). 

In common. Figs. 4-6 reveal a general trend. That is the decrease of the total exchange 
integral J as the inter-rung charge-charge correlations within a ladder increase. We believe 
this trend together with less obvious properties (i)-(iv) are of general significance not re- 
stricted to the particular mechanism we used to describe the CO phase transition in the 
ladder system. 

In the present paper we concentrated on a study of the the intra-ladder superexchange in 
NaV205 . To derive a more complete spin model in the Trellis lattice a thorough analysis of 
the inter-ladder magnetic interactions is required. Such an analysis for the charge disodered 
phase in NaV205 was done in Ref.( |ll|). There the authors claimed that due to the almost 
perfect cancellation of the triplet- and singlet-interactions in a-direction the ladders are ef- 
fectively decoupled to give essentially ID magnetic behavior of the entire spin system. In 
support of this conjecture we recall also that in the charge disordered Trellis lattice the inter- 
ladder Heisenberg spin-spin interactions are strongly frustrated. These arguments, however, 
cannot be transfered directly to a zig-zag CO state. The DFT calculationtO of the exchange 
integrals in a zig-zag CO state have shown, for instance, a rather considerable ferromagnetic 
(FM) inter-ladder coupling0. The opposite sign of this coupling is not surprising in view of 
different symmetry of V-O-V bonds connecting neighbouring ladders. It is a a challenging 
problem to reveal the microscopical origins and estimate a value of the resulting FM coupling 
between neighbouring ladders, which can be done in a properly developed superexchange 
theory taking into account also effects of direct V-V exchange. 
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APPENDIX: 



In this Appendix we exploit an extended model including inter-rung Coulomb interac- 
tions. Our aim is to show that these interactions lead to the symmetry breaking term ~ Ae^ 
in Eq.(^). 

We choose the notation Vab for the Coulomb repulsion between the electrons on the 
diagonal V-V bond, i.e. on a pair of nn V-ions belonging to adjacent ladders (Fig.(l)). We 
assume the parameter Vab to be not too large otherwise the instability with respect to the in- 
line charge ordering would be the dominant one. Now the necessary addition to the single 
ladder model, Eq.(|l]), can be presented in a form Hab = ^{ij} Hal''\ where the summation 
is over the pairs {ij} of nn rungs belonging to adjacent ladders. For a particular {ij} pair 
one has 

i^i?'^ = Vab{^L, jRHnf = Vab{^R, 3L)nfn] (Al) 

where the composite index [i^ means the position of V-ion on the left (right) side of the 
i-th rung; the definitions for the electron density operators n^{a = L, R) are the same as in 
the main text. In the pseudospin representation the term Hab takes the following form 




= Hab,l + Hab,2 (A2) 

where 

= ^ 5^ [Vab{tL,jR) - Vab{tR,jL)] (A3) 

jei 

and summation in (^^) runs over the rungs neighbouring to the i-th rung. 

Without distortion in the Trellis lattice all the diagonal bonds are identical, Vab{iL,jR) = 
Vabi^Ryji) = Vab{\ro\). In this case Ae^ = and the inter-ladder coupling is given by the 
bilinear term Hab,2- Assuming a zig-zag charge ordering in the ladders one may check that in 
the mean field approximation this term is averaged to zero due to a special geometry of the 
ladder stacking. It means that on the mean field level the ladders in the non-distorted Trellis 
lattice are disconnected. Below we argue briefly that a symmetry breaking source appears 
if the lattice symmetry is lowered as reported, for instance, in Refs.( ^, |^) for NaV205 . 

Liidecke et al.@ proposed for the low- T lattice superstructure of NaV205 two species 
of ladders, denoted as A and B in Fig.l. In each A- ladder there is a weak transverse, 
along the a-axis, modulation of the rung positions while the lattice structure of the B- 
ladders remain symmetric as in the high-T phase. The doubling of unit cell means, in 
particular, that in the A-B-A'-B .. array the rung's position modulation in the A- and 
A'-ladders are of opposite signs. In the modulated Trellis lattice a shortening/lengthening 
(±5) of the diagonal V-V bonds leads to a special spatial pattern for the intersite Coulomb 
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parameter variation Vab{\fo =F 6\) ~ Vafe(|ro|) =F V\4b5 = V^*'^'^'* (see Fig.l). By using Eq. (|A3|) 
one can easily check that this pattern produces in the A-ladders an effective longitudinal 
field which is staggered along the ladder b-direction, Ae^^ = (— 1)*A£° , with the amplitude 

A£° ~ 2| VKhf^l- Moreover, the field Ae°^ is also staggered along a-direction, which sustains 
the zig-zag CO state with the A- and A'- ladders being ordered in the antiphase manner. 
Due to smallness of the lattice distortion, |5|/|ro| ~ 10~^, we expect that Ae^ <C Vb for the 
actual values of Vb providing the CO phase transition. 

It follows from Eg. (|A3|) that a symmetry breaking field is not present explicitly in the B- 
subsystem, i.e. Ae% = 0, for the spatial pattern in Fig.l. Strictly speaking, a presence of 
non-zero field Ae^ is not necessary since the long-range zig-zag CO in B-ladders, if it is 
present, can be explained due to spontaneous symmetry breaking if A > Ac. 

Note that besides the extracted staggered field term the inter-ladder bilinear interaction, 
Hab,2i is still present. With a zig-zag ordering in A-/or both in A- and B- subsystems this 
term being considered in the MP- approximation is averaged to zero even in the modulated 
Trellis lattice structure. 
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FIG. 1. Schematic representation of the crystal structure of a vanadium layer of NaV205 . 
The vanadium ions are denoted by open circles, the hopping amplitudes (ta, h and 1^) and Coulomb 
repulsion parameters (V^, V?,) are ascribed to the intra-ladder bonds; the on-site Coulomb repulsion 
U is also shown. Arrows indicate the weak modulation (not in scale) of V-ion positions in the low-T 
phase. The single (double) dashed lines are associated with V^*'^'* (V^'"^ ■*) which are the inter- ladder 
Coulomb repulsion parameters modulated in the low-T phase. 
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FIG. 2. The excitation energies in the singlet (s) and the triplet (r) channels as functions of 
Ae. The model parameters used are (in eV): ta = 0.38, h = t^, = 0.085, U = 4,Va = 0.5. 
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FIG. 3. Pscudospin correlation functions {T?T?_^^g = and {TfTf_^-^^)g ~ as functions of 
A^^=2ia/Vft. The order parameter ruz = \ {Tf)g\ = rj/2 calculated with the exact ID Ising exponent 
/3 = 1/8 is shown with a dashed line. 
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FIG. 4. Variation of the exchange integral J*^*^ and J^'^^ with respect to (the other model 
parameters arc the same as in Fig.l and /3 = 1/8). The integrals are calculated either in the 
modified MF approach (dashed lines) or with the effect of charge fluctuations included (solid 
lines) . 
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FIG. 5. The exchange integrals J^*') and J^'^^ as functions of Vb for various values of Va 



21 




FIG. 6. Variation of the total exchange integral J = J^*^ + J^'^^ in the wide range of the 
varying parameter A^^=2ta/Vf, (a) and in the ordered state as function of the order parameter 
(b). The curves (1) and (2) correspond to two sets of the hopping amplitudes with = i;, > and 
td = respectively (see Sec. V). 
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